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We discuss the calculation of evanescent states in quasi-one- dimensional quantum wires in the 
presence of Rashba spin-orbit interaction. We suggest a computational algorithm devised for cases 
in which longitudinal and transverse motions are coupled. The dispersion relations are given for 
some selected cases, illustrating the feasibility of the proposed computational method. As a practical 
appUcation, we discuss the solutions for a wire containing a potential step. 

PACS numbers: 73.63.Nm 



I. INTRODUCTION 

Quantum wires formed in two-dimensional electron 
gases are electron waveguides allowing transmission only 
along one direction. Due to the lateral confinement, at 
typical electron densities a number of energy subbands 
may be occupied, which are referred to as transverse 
modes or channels. Depending on the behavior of the 
electronic states along the wire they are usually classified 
in propagating and evanescent. The latter decay with 
the distance and, therefore, are irrelevant in the asymp- 
totic regions of the quantum wire, where only propagat- 
ing modes can exist. Nevertheless, evanescent modes are 
of paramount importance for a wire with inhomogeneities 
such as potential barriers or wellsi This is because the 
evanescent modes strongly influence the scattering con- 
ditions determining the reflection and transmission coef- 
ficients from a given potential inhomogeneity. This way, 
the evanescent modes determine the amplitudes of the 
propagating modes in the asymptotic regions of a wire 
with a scattering center. 

The wavenumber k gives the dependence with the dis- 
tance X along the wire for each specific mode as exp(iA:x). 
While for propagating modes /c is a real number, for 
evanescent states it must obviously have an imaginary 
part. In the most general case k may also contain a real 
part for evanescent modes. When transverse and longitu- 
dinal motions are decoupled the wavenumbers for prop- 
agating and evanescent modes can be trivially obtained. 
Indeed, if transverse motion energies are quantized by a 
set of /c-independent eigenvalues with n = 1,2,..., 
the propagating mode wavenumbers fulfill 



±^2m{E-en)/n {E > Sn) 



(1) 



where E is the total energy. On the contrary, evanescent 
modes are characterized by 



k = ±iy/2m{en - E)/h {E < Sn) • 



(2) 



The above scenario of Eqs. ([T]) and ([2j) with analyti- 
cally known wavenumbers, purely real and purely imagi- 
nary for propagating and evanescent modes, respectively, 
breaks down when motion along longitudinal and trans- 
verse directions are coupled. In this case, transverse 



quantization is not given by a fixed set of eigenvalues 
Sn but must be solved explicitly for each E and k. This 
makes, in practice, the calculation of evanescent states a 
formidable task. The mathematical difficulty is exacer- 
bated by the fact that for general complex /c's one looses 
standard properties such as the Hermitian condition of 
the eigenvalue equation, required in most matrix diago- 
nalization algorithms. As a matter of fact, given E the 
problem can be mathematically presented as a non-linear 
eigenvalue equation for k. This situation is found, for 
instance, when a magnetic field is coupled with the elec- 
tron's orbital motion^ or in the specific case addressed in 
this work, when a spin-orbit coupling is present. Evanes- 
cent states with spin-orbit interaction have been recently 
considered in Refs. and|3 in the context of plane- wave 
formalisms, adequate for interfaces between flat bottom 
potentials. However, to the best of our knowledge, a 
general method to obtain the evanescent modes in trans- 
mission channels with arbitrary transverse confinements 
is lacking in the literature. 

We stress that the evanescent states do not fulfill 
proper boundary conditions along the channel since the 
fact that they vanish in one direction necessarily implies 
their divergence in the reversed one. These modes are, 
therefore, not physically realizable in the entire chan- 
nel, for X spanning the interval (— oo,+oo). Neverthe- 
less, they are extremely important because many physi- 
cal states do behave as evanescent in restricted domains. 
This unphysical condition of the extended evanescent 
states explains why, mathematically, they are solutions 
of a non Hermitian problem, while physical states always 
originate from Hermitian operators in quantum mechan- 
ics. 

In this work we shall discuss a practical algorithm to 
obtain the evanescent modes when transverse and longi- 
tudinal motions are coupled through a spin-orbit term of 
Rashba type,^ typical of two-dimensional electron nanos- 
tructures. Although in the specific applications we shall 
assume a parabolic transverse potential, the suggested 
approach will be equally valid for arbitrary transverse 
confinements. Numerical results for selected values of the 
Rashba coupling strength will be presented as an illus- 
tration of the method. Interest in the Rashba spin-orbit 
coupling is mostly due to the tunability of the Rashba 
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strength by means of electric gates which opens a possi- 
bility of spin control in nanostructures, with the Datta 
and Das spin transistor as a well known spintronic device 
proposal.^ 

This paper is organized as follows. In Sec. II we present 
the physical system. Section III is devoted to the practi- 
cal algorithm to compute evanescent states and Sec. IV 
shows selected numerical results and a practical applica- 
tion of our approach. Finally, Sec. V draws the conclu- 
sions of the work. 



II. PHYSICAL SYSTEM 

We consider a two dimensional electron gas lying on 
the xy plane, with a parabolic confinement in the y di- 
rection and perfect translational invar iance along x. The 
corresponding Hamiltonian operator Ho reads 



Px ~^ Py 2 2 

no = + -m^oV 



(3) 



Additionally, a spin-orbit coupling of Rashba type, with 
coupling strength a, is also active 



(4) 



The coupling strength a is assumed constant throughout 
the system, including the asymptotic regions in x and y 
directions. The total Hamiltonian is thus H = Ho -\- Hr 
and we are interested in the solutions of Schrodinger's 
equation for a given energy E 



{H-E)^{x,y,r])=0, 



(5) 



where =t, | is labelling the spin double valued variable. 
Any energy is physically acceptable since the wire spec- 
trum will be continuous for the propagating modes and, 
moreover, since we also consider evanescent modes even 
negative energies will yield solutions to Eq. ([5]). 

As mentioned in the introduction, the spinorial wave- 
function is assumed to be separable in the following form 



ikx 



(6) 



It is convenient to express the spin-dependent part in 
terms of the eigenspinors of Xx±{v)^ 



(l){y,r]) = (l)i{y)xx+{v) ^ My)Xx-{v) • 



(7) 



Equation ([5]) can be recast as a matrix equation for 
the amplitudes (/)i,2(y), 



ho 



la- 



2m 



E 



iak 



-iak 
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ho^^a^,^^ 




i{y) 
h{y) 



0, 



(8) 



where we have defined the transverse oscillator operator 



ho = 



h^d^ 

2m dy'^ 



-mu^y^ 



(9) 



Equation ([8|) is the central equation we intend to solve 
in this work. Although it has the formal appearance of a 
linear eigenvalue equation, with eigenvalue E^ in our case 
the energy is given and what is unknown is the wavenum- 
ber k. Notice also that for complex /c's the matrix in Eq. 
dH]) in non Hermitian, which prevents the use of stan- 
dard matrix diagonalization routines. This invalidates 
the computational strategy normally used for propagat- 
ing modes and consisting in: a) preassign a value to k; 
b) diagonalize ([8]); c) find a posteriori what /c's give as 
eigenvalues the energy of interest E. 



III. THE ALGORITHM 

An essential property allowing a strategy to solve Eq. 
(|8]) is that the amplitudes 0i,2(^) must fulfill the bound- 
ary condition that they vanish for ^ ±oo. Of course, 
^1,2 (^) must also be smooth functions of ^, continuous 



and with continuous derivatives for all y^s. These con- 
ditions can be fulfilled at the same time only for some 
specific wavenumbers representing propagating (real /c's) 
and evanescent (complex /c's) modes. 

Our algorithm is based on the use of finite differences, 
discretizing the y axis in N uniformly distributed points 
in the interval [ymin^Vmax]- The derivatives can then 
be computed using n-point formulas,^ i.e., n neighboring 
points are required to compute the derivative at each grid 
point. We transform our problem into a linear system of 
2A^ scalar linear equations, yielding the 2A^ unknowns 
corresponding to and ^2 on the grid. This linear sys- 
tem has a solution fulfilling the boundary conditions at 
ymin and ymax for any /c, real or complex, but only for the 
physically acceptable wavenumbers both functions have 
continuous first derivative at an arbitrarily chosen match- 
ing point ym- Since the derivatives may be discontinu- 
ous at ym it is essential in the calculation of the deriva- 
tives not to invoke points lying on opposite sides of the 
matching border. This is accomplished with the use of 
non-centered formulas for the derivatives at points near 
the matching border. 

At the matching point we do not impose the discretized 
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FIG. 1: (Color online) Wavefunction amplitudes illustrating 
the algorithm. Real and imaginary parts correspond, respec- 
tively, to (r) and (i) superscripts. The total energy is ^ = /icjo 
while th e Rashb a strength is chosen as a = 0.3/ia;o^o, with 
^0 = yJh/mujQ the oscillator length. The position of the 
matching point is given by the dashed line. Upper and lower 
panels correspond to the indicated wavenumbers. 

Eq. (|8|) but, instead, we ensure a solution different from 
the trivial one (vanishing for all y^s) by arbitrarily choos- 
ing that (j)i{ym) = 1- An additional equation is needed 
to have as many equations as unknowns. We impose con- 
tinuity of the first derivative of (j)2 at ym- In summary, 
the resulting linear system is given by: 

a) (27V — 2) equations obtained by discretizing ([8]) for 
y 7^ Vm^ without crossing the matching border for 
the derivatives, 

b) ^l(^m) 1, 

c) {d(l)2/dy)L - {d(j)2/dy)R = 0, where {d(j)2/dy)L,R 
represent 02 derivatives at ym using noncentered 
formulas with the left (L) and right {R) neighbor- 
ing points. 

The suggested algorithm does not yield normalized 
transverse amplitudes but, of course, normalization 
/ + \4^2{y)\'^dy = 1 can be trivially imposed after 

the solution has been obtained. 

The matrix representing the linear system is greatly 
sparse and can be very efficiently solved using standard 
numerical routines^^ The precision of the method in- 
creases when using larger n's for the derivatives. We 
have used up to 11-point formulas but, typically, 5 or 
7 points already provide quite accurate results. As an 
example. Fig. 1 shows the amplitudes (j)i^2{y) for two dif- 
ferent wavenumbers: upper panel shows a discontinuous 
first derivative at the matching point, indicating that the 
chosen k is not physically valid; the situation is different 
in the lower panel, with perfectly smooth functions. 

Finding the physically acceptable wavenumbers is done 
by exploring the complex-/^ plane. In practice, we define 
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FIG. 2: (Color online) Contour lines of J^(k) showing the 
position of the nodes. We have used E — 0.75h(jJo and a — 
0.3/ia;o^o- 

the real function 




and look for the zeros of JT by sweeping Re(/c) and Im(/c), 
the real and imaginary parts of the wave number, within 
a preselected range. 

IV. RESULTS 

A. Symmetry considerations 

Figure 2 shows the contour lines of circling the posi- 
tion of the nodes for selected values of E and a. Having 
determined the approximate locations, it is then a simple 
matter to zoom in and accurately determine the nodes. 
One must be careful, however, not to miss nodes for 
which the corresponding minimum is very narrow. The 
figure illustrates a symmetry in the k plane: if a mode 
has a certain k then wavenumbers obtained by chang- 
ing sign ±Re(/c) and/or ±Im(/c) are also physically valid. 
Similar conditions were discussed in Ref. [2| for a quan- 
tum wire in a perpendicular magnetic field. Evanescent 
modes with wavenumbers having a nonzero real part will 
therefore come in groups of four, such as the modes in 
Fig. [2] with k = ±0.247 ± 1.166z, corresponding to k, /c*, 
— k and — /c*. 

Mode degeneracies can be explained taking into ac- 
count that if a solution to Eq. (|8]) is characterized by a 
given wavenumber and amplitudes {/c,0i,^2} then, by 
simply taking the complex conjugate of Eq. ([8]), we find 
a solution having {/c*,025^i}- On the other hand, our 
starting Hamiltonian Ti is time reversal invariant which, 
for a spin 1/2 system, implies that the solutions must 
appear in degenerate pairs of time reversed states known 
as Kramers doublets. Using the time reversal operator 
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for a spin 1/2 system, 9 = —icjylC^ where /C represents 
complex conjugation, we find the Kramers partner of 
{/c,^i,^2} as { — A:*, summary, with both 

complex conjugation (cc) and time reversal symmetry we 
obtain the fourfold degeneracy 

el -e: (11) 

{-k\r2,-4>i} A {-fc,-0i,^2}. 

It is worth stressing that the degeneracies expressed 
by Eq. (pTj) are valid irrespective of the ^-inversion sym- 
metry of the wire potential. This is a qualitative differ- 
ence with the result for wires in perpendicular magnetic 
field, ^ where it was found that asymmetric wires display 
asymmetric Re(/c) branches, i.e., they do not fulfill the 
symmetry k ^ —k. The difference originates from the 
breaking of time reversal invariance by the external mag- 
netic field. 



B. Mode dispersion diagrams 

The mode dispersion for the same wire of the pre- 
ceding figures is shown in Fig. [3l We only plot posi- 
tive wavenumbers noting that, with the above mentioned 
symmetry, wavenumber signs can be inverted giving a 
four-fold degeneracy of each evanescent mode and two- 
fold for the propagating ones. The dispersion of prop- 
agating modes shows a familiar picture, already inves- 
tigated in detail^ii^iiiii^ii^ii^ Thcrc is a threshold en- 
ergy En^"^ for the activation of the n-th propagating 
mode. Due to the spin-orbit coupling, when E slightly 
exceeds En^"^ , propagating states with nonzero wavenum- 
ber kn^\ belonging to the n-th mode, are allowed. 

Surprisingly, the evanescent modes in Fig. [3] resemble 
those of asymmetric wires in vertical magnetic field dis- 
cussed in Ref. 2, in spite of the present wires being sym- 
metric and having no field. For E slightly below En^^ 
there is an evanescent mode whose lm{k) approaches zero 

as E approaches En^^ from below but, quite remarkably, 
Re{k) remains essentially stuck at the threshold value 

for the propagating mode, kn^\ This is true even for 
energies much smaller than En^\ Notice that, on the 
contrary, for evanescent states Im(/c) rapidly increases 
when E decreases, indicating a faster decay of the evanes- 
cent mode when the energy separates from the threshold 
En^^ . Strictly speaking, the real part also varies with E 
although it can be hardly noticed to the scale of Fig. [3l 
Figure m shows the mode dispersion for a case of strong 
spin-orbit coupling, when a = hcooio. As compared to 
the weak coupling case (Fig. [3]), the mode branch dis- 
tribution now looks more complicated. Nevertheless, it 
remains true that below the threshold of a propagating 
mode there is an evanescent mode whose wavenumber ini- 
tially coincides with kn^\ the wavenumber of the prop- 
agating mode at threshold. If the energy is further de- 



creased, Re{k) does not remain almost constant as in 
Fig. [31 but displays a rather strong variation, first de- 
creasing and then increasing again. It is worth stressing 
that the symmetry discussed above still applies to Fig. [H 
and one can obtain valid wavenumbers by reverting signs 
in either Re(/c) or lm{k). Another conspicuous feature 

of Fig. [His that en\ the branches giving Im(/c), are not 
monotonously increasing with decreasing energy but dis- 
play a marked wiggle at £^ ~ OAhcJo and O.GhcJo for 63*^ 
and respectively. 

The simpler behavior of the weak coupling case 
suggests an analytical expression for the evanescent 
mode wavenumbers. We have indeed found that when 
a/{huJo^o) ^ 1 the evanescent mode wavenumbers can 
be well approximated by 

k = ±kl'^^ ± i^2m{E^n'''^ -E)/h {E < ^ (12) 

in terms of the threshold energy En^^ and wavenumber 
/cn of the n-th propagating mode. Noting now that in 
the considered limit one has kn^-' ^ kn and En^-' ^ Sm 
where we define the Rashba wavenumber kn = am/h^ 
and the transverse oscillator energies = + l/2)hjOQ^ 
we can further simplify Eq. ([T2|) to 

k = ±kR±i^2m{en - E)/h {E < Sn) . (13) 

This is a very appealing result, stating that in the weak 
coupling case the only modification to Eq. ([2]) introduced 
by the Rashba coupling is to add a constant (energy in- 
dependent) real part given by the Rashba wavenumber. 
This result agrees with the findings of Ref. where the 
density oscillations at the interface between a ferromag- 
net and a semiconductor with Rashba interaction were 
found to have an energy independent wavelength X ^ k^^ 
in the limit of low Rashba coupling. 

Equation ([T3|) can be understood noting that in the 
weak spin-orbit limit one can neglect the Rashba inter- 
subband coupling term, proportional to Pycr^. In this 
approximation the Hamiltonian is diagonal in the Gy ba- 
sis and the corresponding Schrodinger's equation for any 
k implies 

h^k^ , h^kRk ^ ^ 

En + ± — -E = 0, 14 

2m m 

whose solutions for E < Sn and neglecting terms in k\ 
coincide with Eq. ([T3|) . 



C. The potential step 

To illustrate the use of the computed evanescent modes 
with a specific example we have obtained the wave func- 
tions in a quantum wire with Rashba interaction and 
containing a potential step VqQ{x)}^ Electrons incident 
from the left {x < 0) impinge on the step border and. 
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FIG. 3: (Color online) Mode dispersion for a — O.ShuJoio. 
Branch labels indicate mode number n as well as evanescent 
or propagating character as and pn, respectively. In the 
case of evanescent modes, superindexes (r) and (i) are used 
to indicate whether the branch gives the real or the imaginary 
part of the wavenumber. 
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FIG. 4: (Color online) Same as Fig. [3]for a — huuoio. 



assuming their energy is not enough to allow for propa- 
gation in the right region (x > 0), only pure evanescent 
modes will be seen for x > 0. On the contrary, the left 
region will contain three types of modes: incident propa- 
gating, reflected propagating and evanescent. The mode 
wavenumbers before and after the step can be obtained 
from the mode dispersion diagrams, Figs.[3]and[U associ- 
ating to each region an energy measured from the poten- 
tial bottom, i.e., the left Ei = E and right Er = E — Vq 



FIG. 5: (Color online) Density and spin magnetization dis- 
tributions in a wire with a potential step at x = of 
Vb = l.lShuJo and having a = KuqIq. The density contour 
numerical values are given in units of Iq'^ while those of mag- 
netization density are in units of /i/(2^o)- Incidence is from 
the left in mode pi shown in Fig.|4l for an energy E — 0.93/ia;o 
slightly below the threshold for propagating modes ps and pA- 



energies. At the step edge the wave functions must match 
adequately, this condition determining the amplitudes of 
the reflected propagating modes for x < and of the 
evanescent waves for all x. The Appendix details the re- 
sulting equations and their practical resolution method. 
Evanescent states are crucial in this problem for without 
them it is not possible to fulfill the matching conditions 
at X = and for any value of the transverse coordinate y 
(see Appendix). Besides, evanescent states must be con- 
sidered in order to satisfy the obvious requirement that 
the reflected flux coincide with the incident one. 

Figure [5] displays the density and spin magnetizations 
obtained after solving the linear system of equations cor- 
responding to the matching conditions. We have as- 
sumed the same intense Rashba coupling a = fic^o^o of 
Fig. m and a left incidence of unit flux from mode pi, 
with energy E = O.dShcoo^ slightly below the propagation 
threshold for modes ps and p^. The step is chosen as 
Vb = l.lShuJo^o^ which implies that the wavenumbers are 
those of Fig. a at energies (Ei.Er) = (0.93,-0.20)^6^0- 
The spatial distribution of density and magnetization due 
to evanescent modes are clearly seen for x > 0, while on 
the left side there are marked interference effects between 
the coexisting modes. Actually, the latter patterns are 
similar to those discussed in Ref. 4 when considering a 
different physical system; namely, an interface between 
a ferromagnet and a semiconductor. Focussing on the 
evanescent mode side, there is a strong injection of spin 
y magnetization m^, clearly due to the fact that the in- 
cident mode pi is mostly polarized along -\-y (it is not 
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completely polarized due to the admixture induced by 
the Rashba intersubband coupling term apyCix). There 
is also an important accumulation of evanescent rrix mag- 
netization, of different signs on the two sides of the wire. 
Had we considered incidence from mode p2 at the same 
energy, the spin magnetizations would be reverted with 
respect to those of Fig. [5l indicating that the average for 
incidence from the two modes pi and p2 does not pro- 
duce any net magnetization, not even locally. On the 
contrary, the density shown in the upper panel of Fig. [5] 
is the same for incidence from any of the two modes. 

The vanishing spin magnetization when adding the 
contributions from the two left incident channels pi and 
P2 can be interpreted, theoretically, as a manifestation 
of the time reversal symmetry which is conserved by the 
Rashba interaction. Indeed, when Er is such that there 
can be no propagation to the right all incident flux is re- 
flected backwards, yielding a solution which is invariant 
by time reversal. 

A remarkable feature in Fig. [5] is that the density takes 
its maximum value for x > 0, the evanescent mode side. 
This is totally unexpected for purely exponentially decay- 
ing modes. However, since the evanescent wavenumbers 
with Rashba coupling are complex, we may flnd a super- 
position of oscillating and exponentially decaying distri- 
butions, as that of Fig.O This behavior is a peculiarity of 
intense Rashba couplings. To further clarify this we plot 
in Fig. [6] the densities for weak and strong values of a. In 
the weak a case the evanescent density is indeed given by 
a pure exponential decay that can be explained taking the 
wavenumbers given by Eq. (p!3|) . Since Re(/c) is constant, 
it amounts to a common phase for all modes and thus 
irrelevant for the density, remaining only the pure expo- 
nentially decaying contributions. The lower panel also 
shows the density in the strict ID limit, when all trans- 
verse motion is neglected and the Hamiltonian reduces 
to Hid = p^/2m — apxCy/h. The excellent agreement 
with the ^/-integrated 2D density clearly proves that in 
the limit of small a the problem becomes effectively one 
dimensional, as should be expected. 

The relevance of the evanescent modes in the strong 
a case manifests itself in Fig. [6] as conspicuous accu- 
mulations of the integrated density (dashed line), slowly 
decaying towards both sides of the edge. The decay is 
quite sensitive to the values of E and Vq or, equiva- 
lently, to the energies {Ei^Er) flxing the left and right 
modes. This is clearly seen in Fig. [71 which displays the 
^-integrated densities for (Ei^Er) = (0.5, — 0.2)ficc;o and 
(0.5, —0.5)hu;o. Towards the left the decay is now faster 
than in the upper panel of Fig. [6l due to the higher val- 
ues of lm{k) for this Ei. The same argument explains 
the faster decay towards the right for Er = —O.bhuJo. In 
general, we flnd that Ei and Er determine the decay rate 
towards the left and right sides rather independently. 

Most interestingly, the evanescent modes are crucial 
even when the energy Er exceeds the propagation thresh- 
old to the right of the step. This is shown in Fig. [51 where 
the total transmission for left incidence in both modes pi 
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FIG. 6: (Color online) Density along the y — cut (thick 
solid) and integrated in f^^dyp{y)/£o (dashed). Up- 
per panel corresponds to the solution given in the preced- 
ing figure upper panel. Lower panel corresponds to a weaker 
Rashba coupling, with parameters a = O.lhuoio, (Ei,Er) — 
(1.45, 0.45) /ia;o, corresponding to an energy E = IA^Huoq and 
step height Vo = ^^o- The thin solid line of the lower panel, 
almost superimposing on the dashed line, corresponds to the 
strict ID problem. 



and p2 is given for the strong coupling case. We have 
flxed El = 0.93fia;o^o and increased Er above the prop- 
agation threshold of modes pi and p2- In practice, this 
corresponds to lowering the step height Vq = Ei — Er 
while keeping flxed the incident energy E = Ei. The tri- 
angles are the result when evanescent modes are totally 
neglected in subsets {(^n^^n} and {(/)n,^n} of the Ap- 
pendix. There are sizeable differences with the transmis- 
sion obtained when including enough evanescent states in 
both subsets (circles). Quite remarkably, the condition of 
flux conservation which, for the present situation, implies 
that transmission plus reflection amount to a value of two 
unit fluxes (corresponding to the incidence from modes 
Pi and P2) is only fulfllled by the complete calculation. 
Neglecting the evanescent modes leads to a violation of 
this balance, as shown by the inset in Fig. [51 Deviations 
are greater near the propagation threshold where, indeed, 
it is expected that evanescent modes are more important. 



V. CONCLUSIONS 

Evanescent states along a quasi-one-dimensional chan- 
nel couple with the states of transverse motion when spin- 
orbit (Rashba) interaction is active. This greatly compli- 
cates the determination of the evanescent state wavenum- 
bers and wavefunctions since they do not obey a standard 
Hermitian eigenvalue problem. We have devised an al- 
gorithm based on the requirement that the transverse 
spinor wavefunctions fulflll boundary conditions, yield- 
ing a set of linear equations that is solvable for any com- 
plex wavenumber. Physically acceptable wavenumbers 
are those having continuous derivative at an arbitrarily 
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FIG. 7: (Color online) 2/- integrated densities f^^dyp{y)/io 
for the shown values of (Ei,Er) in /icjo units for the strong 
coupling limit a = huoio. Left incidence of unit flux from 
modes pi and p2 is assumed. 




0.0 0.4 

E (^) 



in the following way: 

a) at high energies the mode is propagating; 

b) when the energy crosses the propagation threshold 
the mode becomes evanescent, with a continuous 
evolution of Re(/c), while lm{k) suddenly starts to 
grow from zero; 

c) in the limit of weak Rashba coupling Re(/c) for 
evanescent modes is constant and coincides with 
kR, the propagation wavenumber at threshold, 
while lm{k) is given by Eq. ([2]), the result of un- 
coupled transverse and longitudinal motions; 

d) for intense Rashba couplings a complicated evolu- 
tion of Re(/c) and lm{k) with energy is observed. 

A specific example illustrating the relevance of evanes- 
cent states has been solved, namely the case of a potential 
step in a wire with intense spin-orbit coupling. Density 
accumulations at the edge, as well as distributions of spin 
magnetization have been obtained. Evanescent states are 
necessary to fulfill the boundary conditions at the step 
edge and maintain flux conservation. A remarkable fea- 
ture of the nontrivial evanescent Re(/c)'s is that maximal 
densities can be located on the evanescent mode side. 
We have also analyzed the case when transport to the 
right is allowed, finding an important contribution of the 
evanescent states to the transmissions. 

To summarize, the evanescent modes obtained in this 
work allow the investigation of spin and density distribu- 
tions around inhomogeneities and at interfaces in quan- 
tum wires with extended spin-orbit coupling, a condition 
met in many spintronic devices. 



FIG. 8: (Color online) Transmission for left incidence in 
modes pi and p2 on the potential step when a = huoio 
and El = O.dShuo. The step height for each value of Er is 
Vb = Ei — Er. Circles (solid line) are the complete result while 
triangles (long-dashed) totally neglect the evanescent modes. 
Transmission plus reflection is exactly two in the complete 
calculation and is shown in the inset when evanescent modes 
are neglected. The deviation from the exact value shows the 
relevance of the evanescent mode contributions. 

chosen matching point. 

Evanescent modes in Rashba wires are characterized by 
complex wavenumbers, having real and imaginary parts. 
Due to symmetry, one can invert signs of either Re(/c) or 
Im(/c), or both, and still obtain physically valid wavenum- 
bers. When decreasing the energy the wire modes evolve 



APPENDIX 

In this Appendix we detail the matching conditions 
corresponding to the potential step in a quantum wire 
with Rashba interaction discussed in Sect. IV. C. The gen- 
eral wavefunction is given by a superposition of modes, 
each one characterized by a transverse wavefunction and 
wavenumber {^n(^, v)^f^n}' The potential step VbO(x) is 
such that with the given energy E < Vq for left incident 
electrons only evanescent modes can survive in the right 
side. Let us distinguish the following three sets of states: 
{0n,^n}, the propagating modes incident from x < 0; 
{^n,^n}, reflected propagating and evanescent modes, 
the latter vanishing for x — oo [with negative Im(/c)]; 
and {^n, ^n}, evanescent modes vanishing for x +00. 
The most general state then reads 
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1 tn(^nMy,v)e"'"^ ifx>o, ^ ■ ' 

where a^, and are the usual incidence, reflection and transmission amplitudes. 

The matching conditions at x = require continuity of the wave function and its x-derivative, 

n n n 

YcLnkn(t)n{y,r]) ^Y^nkn^n{y,r]) = ^CnknKiViV) ' (A.3) 

n n n 

Equations (|A.2p and (|A.3p must determine the reflection and transmission coefficients, 6^ and c^, in terms of the 
incident ones and they amount to a linear system of two equations for each value of y and r]. Since we have 
in principle an infinite set of y values, one needs also an infinite set of amplitudes bn and Cn in order to have as 
many equations as unknowns. However, at a given energy the number of propagating states is always finite, which 
is illustrating the fact that one needs to include the evanescent modes, also an infinite set, to fulfill Eqs. ()A.2p and 
()A.3p for arbitrary y. In practice we truncate the sums over evanescent states and and project Eq. IA.2I on 
the set {(j)n} and Eq. lA.SI on {^n}- The resulting linear system reads 

E[4>4>]mnK - [4>4>]mnCn = -Y^i^^lmnt^n , (A.4) 

n ^ — — ^ 

n 

Ekn[^$\mnK - kn[^4>]mnCn = -^kn[^4>]mnan , (A.5) 

n. f ^ n. ' ^ 



where we have introduced the following notation for over- 
lap matrices 

[^^]mn = / ^y^m{y^v)4>n{y^v) , (A.6) 

with obvious extensions for , [cjxp] , [(pcj)] , and [00] . 
Including enough evanescent states {0n} and {0n} the 
calculation converges, giving an increasingly better ful- 
fillment of the matching conditions and of the flux con- 
servation. 

The formalism can be trivially extended to consider the 
case when the energy is such that propagating modes af- 



ter the step {x > 0) are allowed. One just needs to add 
the corresponding propagating modes to the set {0n, ^n} 
while all remaining equations of the Appendix are un- 
changed. 
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